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Advances in measurements of jets and collective phenomena in ultrarelativistic heavy ion colli- 
sions have led to further understanding of the properties of the medium created in such collisions. 
Measurements of the correlations between the axes of reconstructed jets and the reaction plane or 
second-order participant plane of the bulk medium (defined as jet V2), as well as the higher-order 
participant planes (jet v n ), provide information on medium-induced parton energy loss. Addition- 
ally, knowledge of jet v n as well as the ability to reconstruct the event plane in the presence of a 
jet are necessary in analyses of jet-triggered particle correlations, which are used to study medium- 
induced jet shape modification. However, the presence of a jet can bias the event plane calculation, 
leading to an overestimation of jet V2- This paper proposes a method for calculating jet V2 (and by 
extension, the higher jet v n harmonics) and the event plane in an unbiased way, using knowledge of 
the azimuthal angle of the jet axis from full jet reconstruction. 

PACS numbers: 25.75.Bh, 25.75.Ld, 12.38.Mh, 21.65.Qr 



I. INTRODUCTION & MOTIVATION 

The strongly coupled medium produced in high-energy 
collisions of large nuclei at the Relativistic Heavy Ion 
Collider (RHIC) pM] and at the Large Hadron Collider 
(LHC) [SHI] exceeds the energy density at which quarks 
and gluons are expected to be deconfined, and is known 
as the quark-gluon plasma (QGP) [5]. The properties of 
the QGP are studied through many diverse approaches, 
including analyses of (a) collective particle behavior and 
(b) medium-induced parton energy loss [5] . 



A. Collective Flow & Geometry 

Most particles produced in heavy ion collisions result 
from the hadronization of the QGP, and are collectively 
correlated with respect to the initial geometry of the col- 
lision in a way that suggests that the QGP exhibits hy- 
drodynamic flow on the parton level |10j . 

The geometry of the colliding nuclei can be described 
by the "reaction plane," the plane defined by the im- 
pact parameter and the beam direction. However, the 
positions of the participating nucleons fluctuate event-to- 
event so that the matter distribution and its symmetries 
cannot be fully characterized by the geometrical overlap 
of the colliding nuclei. For this reason, the concept of 
participant planes defined by the n th -order symmetries 
of the participating nucleons was introduced [11] [12] . 

As the QGP expands, pressure gradients boost the par- 
ticles along the participant planes, and the coordinate- 
space eccentricity is converted into an anisotropy in 
momentum-space |13j . The azimuthal angular (0) dis- 
tribution of the particles can be expanded in Fourier co- 
efficients with respect to the azimuthal angle of the re- 
action plane (\I/kp), as shown in Eqn. ([!]), or any order 



participant plane (^pp,^,), as in Eqn. ([2| [T3] : 
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In semicentral events the participant zone is roughly 
elliptical in coordinate space, because the geometrical 
overlap region of the colliding nuclei is almond-shaped; 
the minor axis of the ellipse is the 2 nd -order participant 
plane, which is closely aligned with the reaction plane. 
Due to the elliptical shape of the interaction region, the 
V2 term is dominant, in all but the most central collisions, 
when measured with respect to either "Prp or *&pp,2- Re- 
cent theoretical [T5T - H7] and experimental [TBI [T§] work 
has shown that fluctuations in the initial state can also 
result in nonzero higher-order v n coefficients, in partic- 
ular U3, when calculated with respect to the n th -order 
participant plane. Mixed harmonics, such as V3 with re- 
spect to Vl/pp^, will not be discussed in this paper. 

Both the reaction plane and the participant planes are 
theoretical concepts, while the experimental observable is 
the n th -order event plane (denoted by ^EP,n)- Depend- 
ing on the method used to the reconstruct the 2 nd -order 
event plane, the result may be more closely related to the 
2 nd -order participant plane or to the reaction plane. For 
example, it is suggested that asymmetries in the distribu- 
tion of spectator neutrons detected at forward pseudora- 
pidities are more sensitive to the reaction plane geometry, 
while the asymmetries of produced particles at midrapid- 
ity are more sensitive to the participant geometry |20j . 

B. Jets & Jet Quenching 

Partons resulting from hard scatterings in the initial 
stages of heavy ion collisions traverse the QGP and frag- 
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ment into tightly correlated "jets" of particles. The pro- 
duction and fragmentation of these partons is well de- 
scribed by perturbative QCD in proton-proton (pp) col- 
lisions [ST], where no QGP is expected to be formed. 
Therefore by measuring the modification of these jets in 
heavy ion collisions one can learn about parton interac- 
tions with the strongly coupled medium. Suppression of 
particles with high transverse momentum (pt) has been 
observed experimentally in measurements of Raa [2"2T 
[25] and in two-particle correlations [26H32] . and is often 
attributed to medium-induced parton energy loss. 

A few analyses have been performed which merge the 
studies of jets and flow, such as high-px V2 O [34] 
and dihadron correlations with respect to the event 
plane [35] [35] ■ While these methods use high-px particles 
as proxies for jets, recent advances in jet reconstruction 
techniques [37] have made it possible to study highly en- 
ergetic jets in heavy ion collisions directly. 



ysis of jet-triggered particle correlations with respect to 
the event plane. In order to perform such an analysis, 
it is necessary to reconstruct the event plane accurately 
in the presence of a jet. Additionally, for such correla- 
tion studies, knowledge of jet v n is crucial for background 
subtraction purposes. 



II. JET - EVENT PLANE BIAS 

Unfortunately, calculating jet v n is not as simple as 
reconstructingthe event plane using the standard tech- 
nique (Eqn. Q) and then calculating w|j Ct {EP} using 
Eqn. (j5j: 

^ P , n = 1 -t^( ^ WiS ^\ ) (4) 
n Videos [p4i)J 



C. Jet v n 

It is expected that partons produced in hard scatter- 
ings are emitted isotropically in the plane transverse to 
the beam direction, independent of the initial QGP ge- 
ometry. However, it is likely that medium-induced par- 
ton energy loss depends on the length of the parton's 
path through the QGP, and should depend on the QGP 
geometry. Pathlength-dcpcndcnt jet suppression can give 
rise to a difference in the number of jets reconstructed 
in the direction of the event plane and out of the event 
plane, depending on the jet pt range and the parame- 
ters used in jet reconstruction, such as px cuts and the 
resolution parameter R. This effect would result in a 
correlation between reconstructed jets and the n th -order 
participant planes, leading to v n coefficients in a Fourier 
series representation of the azimuthal distribution of the 
reconstructed jet axes (vfjet): 
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Jet v n (denoted by i^ 6 ') describes the correlation be- 
tween the jet axis and the reaction plane or participant 
planes, and should not be confused with a measure of the 
correlation of jet fragments. Non-zero jet v n does not nec- 
essarily indicate that jets experience hydrodynamic flow, 
since pathlength-dependent energy loss is a more likely 
physical explanation. "Jet v n " which describes a corre- 
lation between the jet and the collision geometry, should 
be distinguished from pressure-driven "jet flow," which 
imposes a physical interpretation on any such correlation. 

Pathlength-dependent jet modification could also man- 
ifest itself as a modification of the jet shape (associated 
particle yield and/or width of the jet peaks) as a function 
of the relative angle between the jet axis and the event 
plane. This effect could potentially be seen in an anal- 
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(In Eqn. Q the index i runs over all particles and the 
weights Wi are chosen to maximize the event plane reso- 
lution.) 

If jet fragments are included in the event plane cal- 
culation it will cause the reconstructed event plane to 
be biased towards the jet axis. While the bias may be 
small on an event-by-event basis, it can cause a signifi- 
cant overestimation of the jet v n . Therefore, a method 
for removing jet particles from the event plane calculation 
is desirable. Several ideas for removing the jet particles 
from the event plane reconstruction exist, such as remov- 
ing all the tracks within a cone around the jet axis from 
the calculation, or introducing a pseudorapidity (77) gap 
between the jet and the particles used to reconstruct the 
event plane. The former idea cannot remove all parti- 
cles associated with a jet for reasonable cone sizes, and 
introduces further correlations that contribute to the ar- 
tificial jet - event plane bias. The latter idea requires 
reconstructing the event plane at large pseudorapiditics 
which not only decreases the event plane resolution, but 
also requires a detector with good angular resolution at 
forward angles. 

In this paper a method is presented for calculating jet 
i>2 and the 2 nd -order event plane by accounting for the 
presence of a high-px jet, instead of attempting to remove 
the jet particles from the event plane reconstruction. The 
method is then extended to calculate i4 et and $EP,n- 



III. THE STANDARD EVENT PLANE 
CALCULATION 

The standard method of calculating the event 
plane [38] is by defining a flow vector, Q 2 , that describes 
the bulk particle distribution. The components of Q 2 
are shown in ([6]), where the index i again runs over all 
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the particles used in the event plane reconstruction. The 
vector Q2 and the beam axis define the event plane. 

Q2,x = ^ Wi C0S ( 2< M = cos(2* E p,2) (6) 

i 

Q 2>v = sin(20i) = Q 2 sin(2* E p i2 ) 

i 

Events containing jets, however, can be decomposed into 
two vectors: one describing the bulk distribution (Q2) 
and one describing the jet constituents (A 2 ). When a 
jet is present, the standard event plane method does not 
find the angle of Q2 but rather the angle of the sum of 
these two vectors (G2 = Q2 + A 2 ), as shown in Q: 



G 
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= WiCos(2(f)i) + WiCOs(2(f)i) 

zGbulk iGjet 

= Q 2 cos(2* E p,2) + A 2 cos(2* jot ) 
G2, v ,\ah = Wi sin(20j) 

i 

= ^2 Wi sin(2<^) + 2J Wi sin(2(j>i) 

= Q 2 sin(2* E p,2) + A 2 sin(2* jot ) 

The difficulty is to untangle A2 and Q2 in order to ex- 
tract the angle of Q2 alone. Assuming that jets can be 
accurately reconstructed in heavy ion collisions, there is 
another piece of information that is known: the angle of 
A2, which is the jet axis. 

IV. A NEW EVENT PLANE METHOD 



Solving for u 2 ct and ^ep.2 requires calculating the 



higher moments of G 2 listed in (10 1. Note that terms 
from higher mixed harmonics (on the order of and 
above) are neglected in |l0|): 
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In |9| the quantity Q 2 cos(2(*ep,2 - *pp,2)) is Q 2 ,x 
in the participant plane frame. It is necessary to as- 
sume functional forms for the distributions of <32,ie,pp and 
<52,i/,pp in order to solve the equations in Q and (10 1. 
The distributions of Q 2/x j>p and <32,j,,pp are taken tobc 
Gaussian with standard deviation a . The distribution of 
Q2,x,¥¥ is centered at /i while the distribution of Q2,y,pp 
is centered at zero. The relevant moments of Q2 are 
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The system of equations in ([£])-( 11 ) can be solved to yield 
the parameters fi and a: 
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It is straightforward to solve for jet v 2 using Eqns. dsj|)- 



If G 2<x and G 2tV are calculated with respect to the jet 
axis (instead of in the laboratory frame, as shown in (17])), 
then A 2 only appears in one of the terms: 

G 2 ,xjct = ^2 Wi cos ( 2 (^ - *jot)) (8) 

i 

= Q 2 C0S(2(* E P,2 - *jct)) + A 2 

G 2 ,y,jet = X! Wl sm ( 2 (<fo ~ *jet)) 

i 

= Q 2 sin(2(* E p,2 - * jot )) 

Taking the averages of G 2tX ^ c t and G^yjct over many 
events yields: 

(G^xjet) = {Ql C0S(2(* E R2 - * PP^)))^ + (M) (9) 
(G 2]y j ct ) = 



In (|9j) v 2 ct has replaced the quantity (cos(2(4 r j Ct — 
^pp 2)))- It is assumed in this decomposition that w 2 et is 
independent of Q 2 . 



A. Jet v 2 



The formula for u 2 c in this new method (denoted by 



w J 2 Gt {Q A }) is g iven in E q n - (|13 ) 



jctr n41 _ 4 (G2,a,jet)(G 2)2/ j et .) - 4(G 2 ,x,jetG 2 ^jet) 

(13) 

Note that this method already accounts for the event 
plane resolution by calculating (cos(2(^j e t — ^pp^))) in- 
stead of (cos(2(4 , j Ct — ^ep.2)))- It is not necessary to 
divide the resulting i> 2 et {QA} by the event plane resolu- 
tion. 



Event Plane 



Additionally, Eqns. ([£])-( 12 1 can be solved for (A 2 ), 
the average shift in G 2iX ^ ct due to the jet. 



(A 2 ) = (G 2 ,x,jet) - A*«2 
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(14) 



4 



Once (A 2 ) is known, it can be subtracted from C?2,x,jet 
event by event. Then the event plane can be calculated 
in a way that is, on average, not biased towards the jet, 
using Eqn. (15) which is based on Eqn. Q: 
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£i«i cos (2(&-¥ jot ))- (A 2 ) 



jet 
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data. The thermal (T — 0.291 GeV, which corresponds 
to approximately ^/snn = 200 GeV) background event 
is created such that the 2 nd -order participant plane is 
correlated with \I/j e t according to Eqn. ([3| in order to 
produce a jet v 2 . The pr- and centrality-dependent v 2 
values of the background particles and the multiplicities 
of the background event are chosen to be consistent with 
observations from STAR \MM at RHIC. 



V. HIGHER HARMONICS d, et 



A. Jet V2 



The quantities G„, x ,jet and G n . y j ct 
at any order n: 



can be evaluated 
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However, beyond the definitions of G 
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and G 



the method described above docs not depend on the value 
of n. Therefore the equations listed in (|9|-(13) can be 
used to solve for any tri?* . Solving for the rr h -order event 
plane utilizes Eqns. (14) and (17), and is analogous to 
Eqn. JIB): 
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Jet v 2 is calculated from all charged tracks in the com- 
bined event, without pp-weighting. The results are sim- 
ilar when px-weighting or a track px cut is used. For 
each bin in iV c h, the number of charged particles at mid- 
rapidity in the event, jet v 2 is calculated using the stan- 
dard method (Eqns. j4|- (|5|) and with the new method 
proposed here (Eqn. (|l3"|)). The results of the simulation 



are shown in Fig. 1 for 



Jet _ 



0.0 and 



J et 



0.3 in two 



"2 — u 2 
cases: E^ = 10 GeV and = 30 GeV. 

The results of the simulation indicate that the standard 
method of calculating v J 2 et {EP} leads to an overestima- 
tion of the true jet v 2 by 5-20% (or more, in peripheral 
collisions) due to the jet biasing the event plane calcula- 
tion. However, the new method of calculating u 2 ct {QA} 
obtains accurate values of jet v 2 . The observed discrep- 
ancy between the simulated v 2 et and t>2 6t {EP} is largest 
in the most central events where the bulk v 2 is low and 
the event plane is not well defined, and in the most pe- 
ripheral events where the ratio of jet fragments to bulk 
particles is significant. The artificial jet - event plane 
bias increases with jet energy. 



VI. SIMULATION 

This method has been tested on simple simulated 
events, each of which consists of a PYTHIA [3^1 HO] (PP~ 
like) jet embedded in a thermal background. PYTHIA 
jets are simulated in an interval of ±2 GeV around a 
given transverse energy . The charged tracks in the 
PYTHIA event are analyzed with the anti-fcx algorithm 
from the FastJet package [H] (with a resolution param- 
eter R — 0.4) to obtain the charged jet pt (Px'' ' 1 ) an d 
the angles of the jet axis (\I/jet,??jot)- The charged jet pr 
from FastJet is required to satisfy (X et - 2 GeV) < 



^jet,ch ^ ^_ fgjjeb _|_ 2 GeV^j ; this requirement ensures 

that at least one of the jets from PYTHIA is inside the 
detector acceptance. Only tracks at midrapidity (|ry| < 1) 
are used in the simulation to mimic the acceptance of the 
STAR and ALICE TPCs gHH], and the jet axis is re- 
stricted to lie within |?7j e t| < 1 — R. Note that the jet 
axis and px are reconstructed prior to embedding the 
jet in the thermal background, to avoid issues related to 
full jet reconstruction in a heavy ion environment, al- 
though such a separation is not possible in experimental 



B. Event Plane 

The event plane was calculated using the standard 
method (Eqn. Q) and the new method (Eqn. |l5"| )) and 
then compared with the simulated participant plane in 
Fig. [2] Both methods accurately reconstruct the partici- 
pant plane, and the results are shown separately for the 
cases in which the jet is in-plane (l^jct — ^pp^I < vr/6) 
and out-of-plane (l^jet — ^pp^I > 7r/3). It is clear from 
Fig. [2] that the standard event plane resolution depends 
on the orientation of the jet with respect to the partici- 
pant plane. When the jet is aligned with the participant 
plane, the resolution of the standard event plane method 
is high because the jet pulls the event plane towards the 
participant plane. When the jet is perpendicular to the 
participant plane, the standard event plane is pulled away 
from the true participant plane, thus lowering the reso- 
lution. However, the resolution of the event plane cal- 
culated with the new (QA) method does not depend on 
the orientation of the jet to the participant plane, as is 
desirable. 

Table [I] shows the standard deviation (a) of the event 
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FIG. 1. Jet V2 is calculated using the standard (EP) method (open symbols) and the new (QA) method (closed symbols) in 
simulated events consisting of a PYTHIA jet embedded in a thermal background (T = 0.291 GeV). The simulated jet V2 is 
denoted by a solid line for two cases: v 2 e = 0.0 (left column) and v 1 ^ = 0.3 (right column). The results are shown for two 
jet energies: = 10 GeV (top row) and = 30 GeV (bottom row). The pr- and centrality-dependent V2 values and 

multiplicities of the background are consistent with observations from STAR [44, 45 . Statistical errors are drawn but are often 
smaller than the symbol size. The "jumps" in the distribution of v^jEP} correspond to the edges of the centrality bins in [44] 
that were used to simulate the bulk vi- 



plane from the participant plane, for both the standard 
method and the new (QA) method, in two centrality bins 
and for three orientations of the jet to the participant 
plane. The results show that while the standard event 
plane resolution depends on the orientation of the jet, 
the QA event plane resolution is constant. Furthermore, 
the resolution of the QA event plane does not differ sig- 
nificantly from the resolution of the standard event plane 
in either of the centrality classes. 

C. Jet «3 

The simulation was modified to include a background 
modulated by i>2 and v 3 and a jet correlated to ^pp^ 
and "fpp^ via u 2 et and v 3 ^ 1 . The underlying event i>3 was 
chosen to equal the t>2 measured by STAR in the 5-10% 
centrality bin 44J. The 2 nd - and 3 rd -order participant 
planes are each correlated with the jet axis, but are not 
explicitly correlated to each other. Figure [3] shows the 
results of a simulation in which 30 GeV jets are chosen 
to have =0.1 and i^ ' =0.3. The new jet v n method 
can accurately measure the jet V2 and V3 harmonics in- 



dependently. The standard method estimation of jet v% 
is less biased than the estimation of u 2 et , because in the 
case of odd harmonics, the two jets in a dijet pair pull 
the event plane in opposite directions, and thus the bias 
is largely canceled. 



D. Jets at the LHC 

The simulation was also modified so that the «2 values 
of the background particles and the multiplicities of the 
thermal (T — 0.350 GeV) background event are consis- 
tent with observations from ALICE gHl SZ] at the LHC. 
The results are shown in Fig. for = 100 GeV and 
v 2 ct = 0.3. Due to the higher multiplicities at LHC en- 
ergies, the standard event plane method does not over- 
estimate jet V2 as drastically as at RHIC energies. The 
new method also successfully recovers the correct jet i»2 
at LHC energies. 



G 




FIG. 2. The event plane is calculated using the standard (EP) and new (QA) methods in simulated events consisting of a 
PYTHIA jet with = 30 GeV embedded in a thermal background (T = 0.291 GeV). The p-r- and centrality-dependent vi 
values and multiplicities of the background are consistent with observations from STAR [441 145] . and =0.3. The difference 
between the reconstructed event planes and the participant plane are folded into the interval (— f , f) and are shown for two 
cases: (a) the jet is aligned with the participant plane and (b) the jet is perpendicular to the participant plane. The results 
are integrated over A r c h. 



Jet - Participant Plane Orientation <tep ctqa 


< |*jot - *PP,2| 
7T/6 < |*j ot - fpp,2 
7t/3 < *jct ~ *PP,2 


< tt/6 0.3419 ± 0.0004 0.4378 ± 0.0004 

< tt/3 0.4461 ± 0.0004 0.4374 ± 0.0004 

< tt/2 0.5755 ± 0.0006 0.4397 ± 0.0006 




Centrality 


OEP °"QA 


Peripheral (107 < A ch < 800) 0.3426 ± 0.0004 0.3399 ± 0.0004 
Central (800 < A ch < 1546) 0.4912 ± 0.0004 0.5157 ± 0.0003 



TABLE I. Standard deviation of the event plane from the participant plane, when the event plane is calculated with the 
standard method (oep) and with the new method (<tqa). The results are shown for two centrality bins, and three orientations 
of the jet to the participant plane. Errors are statistical only. 



VII. DISCUSSION & CONCLUSIONS 

As shown by the above simulation study, this new 
method can accurately calculate jet v<i and the 2 nd -order 
event plane in a simple picture of a jet embedded in 
a heavy ion background exhibiting elliptic flow. The 
method also works for higher-order anisotropics. 

This method is tested in a simulation which does not 
include medium-induced jet modification. The decom- 
position in Eqn. (|7|) implies that the bulk particle distri- 
bution is independent of the jet. However, if the jet is 
modified by interactions with the medium, then the bulk 
must be modified to some degree as well. It is expected 
that jet-medium interactions would produce structures in 
the azimuthal distribution of particles that are correlated 



with (or symmetric about) the jet axis, and therefore jet 
quenching would be reflected in the event-by-event mag- 
nitude of A ra , rather than in Q„, to first order. 

This method does not assume anything about the mag- 
nitude of the jet vector (A n ) or about the v n compo- 
nents of the underlying event. However, the equations 
for calculating u^ ot and Vl / EP,n rely on the assumption 
that the fluctuations of Q n<x and Q n ,y in the participant 
plane frame can be described by Gaussian distributions. 
This assumption follows from Ref. [48], although it is 
noted that there are small deviations from the Gaussian 
ansatz in peripheral collisions. However, Ref. [3H] dis- 
putes this assumption, stating that a Gaussian distribu- 
tion is a poor description of event-by-event fluctuations 
in a Monte Carlo Glauber model for all but the most 
central collisions. This new jet v n method may be gen- 
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FIG. 3. Jet V2 and t>3 are calculated simultaneously using the standard (EP) method (open symbols) and the new (QA) method 
(closed symbols) in simulated events consisting of a PYTHIA jet with E^ = 30 GeV embedded in a thermal background (T = 
0.291 GeV). The pr- and centrality-dependent V2 values and multiplicities of the background are consistent with observations 
from STAR 44, 45 . Solid lines denote the simulated (a) jet V2 = 0.1 and (b) jet v$ — 0.3. Statistical errors are drawn. 
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FIG. 4. Jet V2 is calculated using the standard (EP) method (open symbols) and the new (QA) method (closed symbols) in 
simulated events consisting of a PYTHIA jet with E^ = 100 GeV embedded in a thermal background (T = 0.350 GeV) . The pt- 
and centrality-dependent V2 values and multiplicities of the background are consistent with observations from ALICE |46l I47| . 
The simulated jet V2 is 0.3 and is denoted by a solid line. Statistical errors are smaller than the symbol size. 



eralized to account for different distributions of Q n .x.pp 
and Q n .y t pp by rederiving Eqns. (!!)-( 15) with alterna- 



tive non-Gaussian functional forms. 

Although the method presented here assumes perfect 
resolution of the jet axis, accounting for the jet axis reso- 
lution in the measurement of i> 2 et does not require signif- 
icant changes to the equations. This method is improved 
when it is applied in fine bins in centrality (multiplicity) 
and jet energy, in order for the average quantities ((A n ), 
(Qn), etc) to be meaningful. Unfortunately, high statis- 
tics are required in order to obtain reasonable results. 
For this reason the application of this method may not 
be feasible at RHIC, whereas statistics may be sufficient 
at the LHC. 

While jets can influence event plane and V2 measure- 
ments, the production of such high-px jets is rare and 
therefore only a small subset of events are affected in 
most event plane analyses. Furthermore, in u 2 analyses 
of objects that are uncorrelated to jets, the jet - event 



plane bias only contributes to a small decrease in the 
event plane resolution. It is only when calculating the 
V2 of jets themselves, or particles likely to come from jet 
production (such as high-px hadrons), that the jet can 
cause a significant overestimation in the vi calculation. 

A measurement of jet v n for odd n could potentially 
yield more information about jet quenching than the even 
harmonics. Since parton pairs produced in hard scat- 
terings are essentially symmetric under an azimuthal ro- 
tation by 7T, and odd harmonics are antisymmetric under 
the same rotation, it is not possible for jet production to 
have any intrinsic correlation with an odd event plane. 
Barring detector acceptance effects, jet reconstruction al- 
gorithms will likely find (or assign a higher energy to) the 
jet in a dijet pair which undergoes less modification. If a 
nonzero odd w^ et term is measured, it could be indicative 
of the correlation between the n th -order event plane and 
jets which undergo less modification than their recoil- 
ing partners, illustrating pathlength-dependent medium- 



induced jet modification. 



In conclusion, the calculation of jet v n can be beneficial 
for studies of medium-induced jet modification. However, 
reconstructing the event plane in the presence of a jet is 
non-trivial due to the effects of jet fragments on the event 
plane calculation. The artificial jet - event plane bias 
that arises when jet particles are included in the event 
plane calculation can lead to a significant overestimation 
of jet V2- The method presented here utilizes knowledge 
of the jet axis from full jet reconstruction to accurately 
calculate jet t>2 and V3 as well as an event plane that is, 
on average, unbiased by the presence of a jet. 
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